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Abstract 

O . 

We report femtosecond time-resolved reflectivity measurements of coherent phonons in tellurium 
performed over a wide range of temperatures (3K to 296K) and pump laser intensities. A totally 
symmetric Ai coherent phonon at 3.6 THz responsible for the oscillations in the reflectivity data 
is observed to be strongly positively chirped (i.e, phonon time period decreases at longer pump- 
probe delay times) with increasing photoexcited carrier density, more so at lower temperatures. 
We show for the first time that the temperature dependence of the coherent phonon frequency is 
anomalous (i.e, increasing with increasing temperature) at high photoexcited carrier density due to 
electron-phonon interaction. At the highest photoexcited carrier density of ~ 1.4 x 10 21 cm 3 and 
the sample temperature of 3K, the lattice displacement of the coherent phonon mode is estimated 
to be as high as ~ 0.24 A. Numerical simulations based on coupled effects of optical absorption and 

o; 

<3\ . carrier diffusion reveal that the diffusion of carriers dominates the non-oscillatory electronic part 

o i 

of the time-resolved reflectivity. Finally, using the pump-probe experiments at low carrier density 
of 6 x 10 18 cm -3 , we separate the phonon anharmonicity to obtain the electron-phonon coupling 
contribution to the phonon frequency and linewidth. 
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I. INTRODUCTION 



The investigation of ultrafast carrier and phonon dynamics in semiconductors is of spe- 
cial interest as their applications in information technology depends on the understanding of 
fundamental processes like momentum and energy relaxation, as well as mechanisms such as 
carrier- carrier scattering, inter-valley and intra-valley scattering, optical phonon scattering, 
and carrier diffusion [1]. When a semiconductor is irradiated with an ultrashort laser pulse 
with photon energy much higher than the band gap, it excites electrons from valence band 
to the conduction band with excess of energy above the bottom of the conduction band. The 
carrier distribution immediately after the photoexcitation is of non-Fermi-Dirac type (non- 
thermal) {2! where the temperature of the excited carriers is still not defined. Subsequently 
the carrier-carrier interactions lead to a formation of Fermi-Dirac distribution with a definite 
electron temperature (also called as thermalization of carriers) from the initial nonthermal 
state and this process depends on the excited carrier density. At moderate carrier density of 
~ 10 18 cm -3 , this carrier-carrier interaction time is ~ 200 fs in GaAs 2J, |3| which can reduce 
to a few fs at higher photoexcited carrier densities. Following this, the hot carriers diffuse in 
the material as well as emit phonons. The magnitude of cooling time again depends on the 
number of excited carriers and the initial temperature of the distribution. Typically, carrier 
cooling times are of the order of a few picoseconds in semiconductors {4]. There is, finally, 
the electron-hole recombination which happens on the time scale of several picoseconds to 
nanoseconds. 

In recent years, the femtosecond pulses have been used to generate and probe coher- 



ent phonons in solids [5|-|9|]. The excitation mechanism of coherent phonons was thought 
to be impulsive stimulated Raman scattering (ISRS) [5| in transparent materials whereas 
displacive excitation of coherent phonons (DECP) is the dominant mechanism in opaque 
. Later, it was shown that DECP is a special case of ISRS when excited reso- 
nantly 9J. In a resonant ISRS/DECP process, the photo-excitation of electrons changes the 
equilibrium positions of the atoms; the atoms then oscillate around their new equilibrium 
positions. In absorbing samples like Te, the Dember field arising from different diffusion 
coefficients of photo-excited electrons and holes can also contribute to the generation of 



coherent phonons 



10j-[l2|. 



Till now, there have been only a very few ultrafast coherent phonon studies at photoex- 
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cited carrier densities (PCD) as high as 10 20 cm 3 13l-ll9|. A red shift of frequency of the 



totally symmetric A\ phonon was observed in a narrow band gap semiconductor, tellurium, 



X 




LO 21 
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~ 3 at room 



14] . This was at- 



temperature in a degenerate pump probe reflectivity experiments 

tributed to the bond weakening by the carriers {^(J. In addition, the phonon time period 
was found to decrease with increase in the delay time between pump and probe pulses, lead- 
ing to asymmetric line-shape in frequency domain [3]. This linear sweep in the frequency 
with the pump-probe time delay can be termed as phonon chirping which originates from 
the rapid change of photoexcited carrier density in the penetration depth of the sample due 
to carrier diffusion and its interaction with the coherent phonons. Similar phonon chirping 
has been seen in degenerate pump probe experiments on bismuth with PCD of > 3 x 



10 cm . Two pump pulse experiments 
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calculations 
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19] along with first principal density functional 



19] have established that the chirping arises from varying carrier density 



along the probed length in the sample and not due to lattice anharmonic effects. 

We note here that all the previous studies on tellurium have been performed at room 
temperature. Here, we present femtosecond time-resolved reflectivity measurements of co- 
herent phonon dynamics in tellurium performed over a wide range of temperatures (3K to 
296K) and carrier density levels. The non-oscillatory part of the data associated with the 
strength of photoexcited carrier density is found to be dominated by diffusion of the carri- 
ers as revealed by our numerical simulations. The oscillatory part shows highly positively 
chirped Ai phonons excited at photoexcited carrier density of ~ 1.4 x 10 21 cm -3 and the 
chirping is found to be largest at lower temperatures where the carrier diffusion time and the 
strength of photoexcited carrier density were found to be maximum. Apart from the phonon 
chirping, the associated lattice displacement is estimated to be as high as ~ 0.24 A at 3K. 
The other important result is the observation of increased softening of A\ mode with the 
pump fluence at 3K (-0.16 THz/mJ-cm -2 ) as compared to 296K (-0.08 THz/mJ-cm~ 2 ) con- 
sistent with the increased photoexcited carriers at lower temperatures. Using time-resolved 
reflectivity experiments performed at low carrier density levels of ~ 6 x 10 18 cm -3 where 
phonon chirping is found to be absent, the effects of phonon anharmonicity and high carrier 
density are separated and analysed. 
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II. SAMPLE DETAILS 



In our experiments, we have used a polycrystalline tellurium crystal (6x4x2 mm 3 ) 
mounted on a continuous helium flow optical cryostat. Tellurium is a narrow band semi- 



conductor with a band gap of ~ 0.33 eV [15j. It crystallizes into P3i21 structure with the 
point group D% and three atom basis arranged around a screw axis parallel to the c axis 
of the hexagonal unit cell forming a helical structure. The equilibrium lattice constants are 
a = 4.456 A (distance between two helices) and c = 5.927 A. The interatomic spacing in 
a chain is 2.834 A and the helical radius is u = 1.173 A. Group theory gives the allowed 



Raman modes 



2l| as Ax (~ 3.6 THz)+ E^ G (~ 2.8 THz) + E 2 G (~ 4.2 THz). The totally 



symmetric Ai mode corresponds to symmetric inter-chain dilation and compression normal 
to the chain axis, involving bond length and bond angle distortions. In this mode, a and c 
do not change whereas only helical radius u changes. 

The linear absorption coefficient j^j, a, of Te at 1.57 eV is 2 xlO 5 cm -1 which gives the 
penetration depth £ (~ 1/a) of pump light to be ~ 50 nm. At maximum pump fluence of 
3.3 mJ/cm 2 , the excited carrier density is Nq = Fa(l — R)/E p is ~ 1.4 xlO 21 cm -3 which is 
0.8 % of the total density of electrons in the valence band. Here, F is the pump fluence in 
J/cm 2 , R is the reflectivity (0.46 at 1.57 eV) of tellurium and E p is the photon energy in J. 



III. EXPERIMENTAL DETAILS 



Degenerate femtosecond pump-probe experiments in the reflection geometry were carried 
out using femtosecond pulses derived from Ti: Sapphire amplifier producing 50 fs pulses 
with the central photon energy of 1.57 eV at a repetition rate of 1 kHz. The pump beam 
was modulated at 393 Hz with a mechanical chopper and the reflected probe intensity was 
recorded with a Si-PIN diode and a lock-in amplifier. The pulse width at the sample point 
was 65 fs (FWHM) measured using a thin beta-barium borate (BBO) crystal. The spot 
sizes (half width at 1/e maximum) of the pump and the probe beams were kept at 600 /xm 
and 400 /im, respectively at the overlap of the two beams at the sample. The difference 
between the probe intensity reflected from the sample and that of a reference beam was 
recorded as a function of delay between the pump and probe pulses. Both the pump and 
probe beams were kept close to the normal incidence with their polarizations perpendicular 
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to each other. This detection scheme can detect only Ai mode for a single crystal. However, 
the non-symmetric E modes can be seen as well for a polycrystalline sample with major 



crystal axes oriented at an angle from the surface [U|, [17| • In our experiments, we observed 
that the damage threshold is ~ 4.5 mJ/cm 2 which is much smaller than the value of 12 
mJ/cm 2 reported earlier [kJ. The pump ffuence was varied between 3.3 mJ/cm 2 and 0.3 
mJ/cm 2 whereas the probe ffuence was kept at 0.04 mJ/cm 2 . The time-resolved reflectivity 
was recorded at four temperatures - 3K, 80K, 160K and 296K. At the lowest pump ffuence 
of 14 /iJ/cm 2 (probe ffuence ~ 0.8 /ij/cm 2 ), the transient reflectivity data was recorded 
as a function of temperature from 296K to 3K using femtosecond pulses derived from the 
oscillator (50fs, Tsunami, 76MHz, Spectra Physics Inc) in the fast scan scheme j^j]. In this 
scheme the pump delay line is modulated with a shaker at 65 Hz and the time-resolved signal 
was detected with a fast A/D converter (AlXscan data acquisition system, AMO GmbH, 
Germany Inc.) which helps in capturing relative signal changes of ~ 10~ 7 without using a 
lock-in amplifier 231 ]. 



IV. RESULTS AND DISCUSSION 

The summary of results at 3K for various pump fluences is presented in Figfj] where the 
panel (a) is for the time domain data and panel (b) shows the Fast Fourier Transform (FFT) 
of the oscillatory part of the data. The signal contains damped oscillations associated with 
the generated coherent phonons overlaid on an exponentially decaying background. This 
background can be assigned to direct changes in electronic susceptibility by the photoexcited 
carriers j3]. The vertical dashed line in frequency domain is a guide to the eye to appreciate 
the red shift of the peak frequency with increased pump ffuence. The time domain data at 
the lowest pump ffuence of 0.25 mJ/cm 2 can be fitted with the following equation 

Ai? 

— - = y + A el exp(-t/r el ) + B exp(-t/r) cos(2nut + 4>) (1) 
ri 

where r e i is electronic relaxation time; and A e i = (^)ei is the strength of electronic re- 
laxation contributions; B = (^r) p hi t, v, and <fi denote the amplitude, relaxation time, 
frequency, and initial phase of the coherent phonon oscillations, respectively. Here yo de- 
notes the background arising from long relaxation time constant (>1 ns) which may be 
associated with slow electron-hole recombination time. At higher pump intensities, the fit 
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to the data using Eq. <^ was found to be unsatisfactory as the frequency rendering a good 
fit at short time delays deviates considerably from the data at longer time delays. This is 
demonstrated in Fig. [2] (a) and (b) for the time domain reflectivity data and FFT of its 
oscillatory part at 3K and pump fluence of 3.3 mJ/cm 2 where thin grey line is the fit using 
Eq. ([T]). The experimental data is best fitted by taking into account of phonon chirping. 
The transient reflectivity data is modelled as 

— = y + A el exp(-i/r e ;) + ^ B i ex P(-*/ r i) cos(2m/it + Pit 2 + (2) 

8=1,2,3 

where /3 is the chirp parameter describing linear sweep in the phonon frequency with the 
pump-probe delay time. Here the fit is carried out with three modes (represented by sub- 
script, % — 1, 2 and 3) guided by the FFT (Fig. |2](b)). The coherent phonon modes observed 
at 2.45 THz, 3.25 THz and 4.13 THz are attributed to E^ G , h\ and E^ G , respectively. It 
is found that the Eq. (T5]) furnishes a very good fit (black line in Fig. |2^a) and [2]^b)) to 
the data. The chirping is significant for only the Ai mode {% — 2) for which the obtained 
fitting parameters are /?2 = 0.55 ps -2 , z/2 = 3.24 THz, T2 = 0.81 ps and B2 = 0.22. The 
value of Ai phonon oscillation amplitude, B2 is 34% of the maximum change in reflectivity, 
(^f)max( = 0.648). We note that the non- symmetric E- modes were observed only at 3K 
using pump-fluence of 3.3 mJ/cm 2 . For pump fluences lesser than 3.3 mJ/cm 2 , data are 
fitted using Eq. fl2]) by taking only one coherent phonon (Ai). The fits as shown in time 
and frequency domain in Fig. [T]are excellent. 



The lattice displacement of 
materials like bismuth from 16 



;he coherent phonon modes can be estimated for absorbing 



24 1 

(3) 



tt2 3.8 x 1(T 3 Bi F 



gui\e\ 



2e 2 ' 



where Uj is in Angstrom (A), F is the pump fluence in mJ/cm 2 , g is the density of the 
material in amu/A 3 and e is the dielectric constant (e = S\ + je^, ^ P h is the energy of the 
phonon in eV and D = with E as the photon energy in eV. Now, for Te at 1.57 eV, 
D = iff ~ lO-^y- 1 , fjg^ ~ l0 3 eV-\ and e = 32 + jll = 34) at 300K |2^. The 
value of \e\ at 10K is 27 [22]. Since the exact temperature dependence of \e\ is not available, 
we have used |e| = 34 at all temperatures. This will introduce a maximum error of ~ 10% 
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at 10K. The final expression for the lattice displacement in this case is 



QVi\e\ 

This gives a very high value of U2 ~ 0.24 A which is almost twice of ~ 0.13 A estimated in 



Bi 



17]. In terms of the lattice constant a = 4.4561 A, U2 ~ 0.05a which is unusually very 



high. The oscillation amplitude of E^ and E^ are, B x = 0.03 and B 3 = 0.01 respectively 
corresponding to lattice displacements of Ui ~ 0.11 A and U3 ~ 0.05 A. Thus the huge 
lattice displacements caused by the high carrier density excitation may be responsible for 
observing the forbidden Eto modes in isotropic reflectivity detection scheme 17J]. The 



possibility of these E modes due to major crystal orientation in the polycrystalline sample 
is perhaps negligible as these modes were seen only at 3K with the highest pump-fluence. 
Since E-modes are absent at all lower pump fluences and temperatures higher than 3K, 
we proceed now to discuss only the Ai mode's behavior. From here, we make a change in 
notation for the fit parameters corresponding to phonons as, B p h, r p h, v v h, (3 and (p. The 
fit parameters of the non-oscillatory electronic part and the oscillatory phononic part are 
presented in the following sections. 



A. Electronic part 

Fig. |2](c) and [21(d) show the dependence of A e i and T e \ on PCD at four different tempera- 
tures. It is seen from Fig. [2(c) that the strength of the electronic contribution, A e \ increases 
with the pump fluence at all temperatures. A e i is fitted (solid line in the figure) with A e i ~ 
N p where p = 1.5 at 3K and 80K; and p = 1.0 at 160K and 296K. The photoexcited carrier 
relaxation time, r e i increases with the PCD (N) (Fig. 121(d)) and saturates at ~ 0.6 x 10 21 
cm -3 . Figs. [21(e) and[21(f) display A ei and T d vs temperature at different PCD (legend C1-C5 
corresponds to 0.1, 0.2, 0.4, 0.6 and 1.4 x 10 21 cm -3 ). Thick lines in Fig. [21 (e) are the 
linear fits (A e /(T) = A e /(0)- bT). It is seen from Fig. [21(e) that A e ; increases with lowering 
of the sample temperature, reaching a maximum value of 0.36 at 3K (at 1.4 x 10 21 cm -3 ), 
almost 4 times higher than its value at 296K. The trend is same at other PCDs. We do not 
yet understand the temperature and PCD dependences of A e /. The dependence of r e i (Fig. 
12(f)) on temperature is rather anomalous: the electron relaxation time is seen to increase 
with increase in temperature. The solid lines in Fig. [21 (d) and (f ) are drawn as guide to the 
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eye. 

We now try to understand carrier dynamics using the diffusion of photoexcited carriers 
with time. When a femtosecond pump pulse irradiates the sample, it creates a carrier 
density, N on the surface depending on the absorption coefficient and the fluence. The 
carrier diffusion leads to a rapid decay of these carriers across the surface. The diffusion 
coefficients of photo-generated electrons and holes can be different leading to a Dember field 
with a build up time of 100 fs to 200 fs after the pump excitation. Since the carrier 



relaxation times in our experiments are long (2 ps to 6 ps as shown in Fig. EJ^d)), we work 
with the assumption of ambipolar diffusion of carriers, described by the diffusion equation 



26] 



^4 t ^*V-« (5) 

where t p is the pulse width, D a is the ambipolar diffusion coefficient and z is the depth 
inside the medium. We solved this partial differential equation numerically for the carrier 
density N(z,t) with the initial condition N(z, 0)= iVo and the boundary conditions, D a ^.=0 
at z = and and N(L,t) = Nq e~ L ^. The phonon displacement, x p h, given by a damped 
harmonic oscillator where the displacement and the angular frequency, Qq(N) depend on 
the carrier density N(z,t) [l5] 

^ + n (NY[x ph - x (N)\ + = o (6) 

Here T is a damping constant for the phonon. In linear approximation (which is expected to 
be valid since the PCD is less than 1% of the valence electron density), the carrier density 
dependence of phonon displacement, xq(N) and the angular frequency, f2o(iV) are taken as, 
x (N)= x eq + (dx /8N)N and £Iq(N) = Vl eq + (dil /dN)N; x eq and Vl eq are the lattice 
displacement and angular phonon frequency at equilibrium (at low carrier density levels). 
We have used x eq = 0.2686 a (ref. 15j) and the experimentally relevent values {Q eq = 23 
THz and Y = 0.56 THz and [f^] = 1.74 x 10~ 21 THz-cm 3 at 160K [see inset of Fig. 1 
(d)]). Eq. was solved numerically using fourth order Runge-Kutta method with the 
integration step of 2 fs. The dielectric constant e — E\ + je 2 depends on x p h and N(z,t). 
Since e\ = 32, much larger than e 2 (= 11), the changes in e± will dominate the change in 



S 



the reflectivity. Thus, 



AR 



(7) 



R 



? — 



with 




(8) 



where ^ and ? are taken as adjustable parameters to fit the experimentally observed 
background and initial oscillations in the reflectivity As an example, the data at 160K 
and the fit using the Eqs. (El)-® are shown in Fig. [3] in panel (A) for different values of 
pump fluence. The variation of ambipolar diffusivity, D a extracted from the fit are used to 
calculate the diffusion time ~ (, 2 /D a . The panel (B) of Fig. [3] displays thus derived 
(shown by stars) r d vs N (Fig. Bj (a2)-(d2)) and r d vs T (Fig. S (e2)-(h2)). Thick lines 
are drawn as a guide to the eye and they are not due to digital smoothing. Comparison of 
Fig. [3] (a2)-(d2) with Fig. |2] (d) shows that the values of T d and T e \ are comparable and 
their dependence on N is similar. Similarly, it is seen that the behavior of t& vs T (Fig. [3] 
(e2)-(h2)) is compatible with r e vs T (Fig. 12(f)) except at the lowest pump fluence. From 
the above analysis, it is clear that the carrier diffusion plays a major role in determining the 
reflectivity changes in tellurium after the the pump pulse excitation. 

B. Phononic part 

The dependence of phonon amplitude (B ph ) on pump fluence (Fig. 0] (a)) is similar to 
that of A ei (Fig [2] (c)). A plot of B ph vs A d given in the inset of Fig. 0] (a) shows a linear 
dependence between them, implying that the photoexcited carriers contribute mainly to 
the generation of coherent phonons as expected in a DECP process. The phonon lifetime, 
r p h decreases with increasing the pump fluence (see Fig. @] (b)) due to electron-phonon 
interaction. Thick lines in (a) and (b) are the guide to the eye. At 3K, the phonon lifetime, 
r p h is ~ 5.73 ps at the lowest PCD of ~ 0.1 x 10 21 cm~ 3 (corresponding to the phonon 
linewidth 7 = — !— = 0.056 THz); this decreases by 87 % to 0.74 ps at the highest PCD of 
1.4 x 10 21 cm -3 . A similar behavior is seen at all other temperatures. 

Next, we turn our attention to the chirp parameter (3 shown in Fig. H] (c). The chirp 
parameter is fitted (solid line) with ~ JV with p = 2.2 (3K), p ~ 2.5 (80K and 160K) 
and p = 3.8 (296K). The phonon chirping increases with increasing pump fluence; since 
A e i ~ N p , f3 ~ Agj , p . It can be seen from Fig. H] (c) that (3 increases on lowering the 
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temperature. For example, at highest PCD of 1.4 x 10 21 cm 3 , (3 ~ 0.55 at 3K as compared 
to ~ 0.2 at 296K. 

We now look at the dependence of the observed phonon frequency as a function of photo 
carrier density at two temperatures, 3K and 296K as shown in FigJH(d). Here, the variation 
of the frequency is linear with N (linear fits are shown by solid lines in Fig. [5] (d)), as 
taken before in solving the Eq. ([6]). ^(THz) = z/o(THz)-y(THz cm 3 ) A(cm~ 3 ) with uq 
= 3.75 THz and 3.64 THz at 3K and 296K respectively; and y = ± [fjja] = 3.73 x 10" 22 
THz-cm 3 at 3K and 1.86 xlO -22 THz-cm 3 at 296K. The slope y determined here corresponds 
to the change in frequency with pump fluence as ~ -0.16 THz/mJ-cm -2 (3K) and ~ -0.08 
THz/mJ-cm -2 (300K). The value at room temperature agrees well with the theoretically 



estimated value 



151 ] of -0.085 THz/mJ-cm 2 at 296K for tellurium. The modulus of the 



slope, | y | from this linear fit is plotted as a function of temperature in the inset of Fig. [4] (d) 
showing that the softening of phonons with carrier density decreases linearly with increasing 
temperature. 

The lattice displacements of Ai mode using Eq. (pfl) at various PCD vs temperatures are 
summarized in Fig. H] (e) where the solid lines are connected through the data points as 
guide to the eye. Consistent with more photoexcited carriers at 3K (pump fluence of 3.3 
mJ/cm 2 ), U/a is ~ 5 %, which gradually decreases to 2.5 % at 296K, i.e a decrease by a 
factor of 2. The trend is similar for other pump fluences (see Fig. 0] (e)). 



C. Low pump fluence experiments 

To separate the effects of electron-phonon and phonon-phonon anharmonic interactions 
in the frequency and lifetime of coherent phonons presented above, we have performed the 
degenerate pump-probe reflectivity experiments with very low pump fluence of 14 /xJ/cm 2 
(Nq ~ 6 x 10 18 cm -3 ). The varying non-oscillatory electronic background in the data is 
removed using FFT digital filter smoothing to extract the coherent phonon oscillations. Thus 
obtained data was fitted using Eq. (JTJ as the chirping is absent. The typical reflectivity data 
after FFT filtering at 3K along with the fit are shown in Fig|5] (a). The fitted parameters- the 
coherent phonon amplitude Bp^, damping parameter j p h (= ^77^) and the phonon frequency 
u p h are plotted as a function of temperature in Fig[5](b), (c) and (d). It has been shown earlier 



27| that the coherent phonon amplitude, B p h behaves quite similar to Raman peak intensity 
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as L ~ [2n(u P fr/2)+i] ' where Ti{y p h) is the BoseEinstein factor. The fit (solid line) using this 
expression is shown along with the data for B p h in Fig. [3(b) where, B p h is normalized with 
respect to its value at 3K. ^ vh and v are fitted (solid line in Fig. [5] (c) and (d)) with the 



following well known functions 



28 



29| | based on cubic anharmonicity where the phonon 



of frequency v decays into two phonons of equal frequency: 7 p h(T) = 7o + A^ anh {T) with 
A lanh (T) = C[l + 2n{u /2)] and u ph (T) = u + Av anh (T) with Au anh {T) = A[l + 2n{v Q /2)\ 
where uq, A, C and 70 are the fitting parameters (A and C are the measures of third order 
cubic anharmonicity). The parameters obtained from fitting are v = 3.70 ± 0.01 THz, A 
= -0.022 ± 0.002 THz, 7o = 0.010 ± 0.003 THz and C = 0.042 ± 0.003 THz. Fig. Mb), (c) 
and (d) show that the fits to the anharmonic contributions are excellent. 



D. Electron-phonon coupling 

In doped semiconductors, the electron-phonon interactions contribute to the linewidth 
(7i) of the i th phonon of frequency and degeneracy gj which is related to the dimensionless 
electron-phonon coupling constant Xep and density of states at the Fermi level (DOS(ep)) by 
the Allen's formula |3oj |: 7$ = (2tx j gi)XEpjyfDOS(Ep). The question arises what happens to 
phonons due to the photoexcited carriers under femtosecond pulse excitation. We address 
this question and extract the contribution of electron-phonon interaction to the phonon 
frequency and linewidth as described below. 

The temperature dependance of u p h and j p h along with their linear fits at various pump 
fluences are drawn in Fig. [6] (a) and (b)(legends C1-C5 correspond to 0.1, 0.2, 0.4, 0.6 and 
1.4 x 10 21 cm" 3 ). We note that the effect of the carrier diffusion is explicitly taken into 
account while extracting the phonon parameters, v p h and 7^ (shown in Fig. [6] (a) and 0(b)) 
by fitting the total ^ with Eq. [3 It can be seen that the normal anharmonic behavior of 
phonons (i.e, frequency increasing with lowering temperatures) at the lowest pump fluence 
of 0.1 x 10 21 cm -3 changes gradually to an anomalous behavior (frequency decreasing with 
lowering temperatures) as the pump fluence is increased, suggesting a strong electron-phonon 
coupling at higher carrier densities. No such anomalous behavior is noticed in the phonon 
damping term (Fig. [6](b)). The phonon frequency v p h(T) and damping j p h{T) thus deduced 
have contributions both from electron-phonon and phonon-phonon anharmonic interactions. 
The temperature variation of frequency and linewidth of a phonon mode can be thus written 
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as [23,1211 



u ph {T) = u + Au anh (T) + Av e i-. ph (T) 

l P h{T) = 70 + Manh{T) + Aj d _ ph (T) 



(9) 
(10) 



where v e \- v h and ^ei-ph are due to electron-phonon interaction. 

For the low-pump fluence case of 14 /i J/cm 2 (data shown in Figs. [5] (c) and (d)), the 
contribution of electron-phonon coupling is negligible and the temperature dependence arises 
predominantly from anharmonic interactions. Thus, by subtracting the values of v p h(T) and 
lph{T) shown in Figs. [5](d) and[5](c) from the corresponding values shown in Figs. [6] (a) and 
M (b), we obtain only the electron-phonon contribution to the phonon parameters (shown 
in Figs. M (c) andO (d)) at different carrier density levels. The solid lines in Fig. O (c) are 
fit to Av e i^ p h(T) = Au e i^ p h(0) + A(N)T, where the slope A is a function of photoexcited 
carrier density N. It can be seen that both Av e i_ ph and A^ d _ vh strongly depend on the 
carrier density N. On the other hand, A7 e ;_ p / l is weakly dependent on T. More theoretical 
understanding is needed to explain the trends seen in Fig. [6] (c) and (d). 

V. CONCLUSIONS 

Time-resolved reflectivity measurements of tellurium have been performed over a wide 
range of temperatures and photoexcited carrier densities. The relaxation time associated 
with the carrier diffusion increases with the pump fluence which is quantitatively understood 
based on a diffusion model. The temperature dependence of carrier relaxation time, r e i, is 
observed to be anomalous. The oscillatory part of the time-resolved differential reflectivity 
shows highly chirped Ai phonons excited at high carrier densities and the chirping is found 
to be largest at the lowest temperature where the carrier diffusion time and the strength 
of photoexcited carrier density are maximum. Apart from the chirping seen in phonons, 
the lattice displacements associated with them is estimated to be very high (e.g ~ 0.24 
A at 3K). Another important result is the observation of increased softening of A\ mode 
with the pump fluence at 3K (-0.16 THz/mJ-cm~ 2 ) compared to 296K (-0.08 THz/mJ- 
cm~ 2 ) consistent with the increased photoexcited carriers at low temperatures. Using low 
pump fluence transient reflectivity experiments, the effects of phonon anharmonicity and 
the contribution of electron-phonon interaction to the phonon frequencies and linewidths 



12 



at high pump fluence were separated. The temperature and photoexcited carrier density 
dependences of electron-phonon contribution to the mode frequency and linewidth needs to 
be understood theoretically. 
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FIGURE CAPTIONS: 

Figure 1: (a) Normalized time-resolved reflectivity data (open circles) at T=3K at various 
pump fluences along with the fits (line) according to Eq. ([2]). (b) Corresponding FFT spectra 
of the oscillatory part of the data and fit. The vertical dashed line in (b) is the guide to the 
eye to appreciate the red shift of the phonon peak frequency. 

Figure 2: (a) Normalized time-resolved reflectivity change (^r) at T=3K at maximum 
pump fluence of 3.3 mJ/cm 2 along with the fit.(b) The FFT spectra of the oscillatory part 
of time domain part. Open circles represent the experimental data and thick line is the fit 
according to Eq. ([2]) with = 0.55 and the grey line is with (3 = 0. (c) Carrier density 
dependence of strength of photoexcited carriers (A e ;) and (d) photoexcited carrier relaxation 
time {r e i) at different temperatures (T=3K (filled squares), 80K (filled circles), 160K (filled 
triangles), 296K (filled inverted triangles)). The fit, A e ; ~ N p is shown in (c) with p 
=1.5 (3K), 1.6 (80K), 1 (160K and 296K). The lines in (d) are guide to the eye. (e) The 
temperature dependence of A e i and (f) r e i at given carrier densities (CI = 0.1 x 10 21 cm -3 , 
C2 = 0.2xl0 21 cnT 3 , C3 = 0.4x 10 21 cm" 3 , C4 = 0.6x 10 21 cm" 3 , C5 = 1.4x 10 21 cm~ 3 ). The 
lines in (e) are the linear fits and in (f) are guide to the eye. 

Figure 3: (A) Time domain data at 160K (open circles) and the corresponding fit (line) 
using the diffusion model (see text) at different fluences (al)-(el). (B) The diffusion time 
(filled stars), r d = £ 2 /D a versus N (a2)-(d2) and T (e2)-(h2) (see text). Here D a is the 
carrier diffusion coefficient and the solid lines are drawn as guide to the eye. 

Figure 4: PCD dependance of phonon fit parameters: (a) phonon amplitude (B p h), (b) 
phonon life time {T p h), (c) phonon chirping parameter (/3) along with a fit ~ N p (2.2 < p > 
3.8), and (d) the phonon frequency, v p h at two temperatures-296K and 3K. The inset of (d) 
shows the dependance of the slope (filled stars), y (= ^rf^y THz cm 3 ) of the linear fits to 
the data as function of temperature (solid line is the linear fit), (e) The variation of the 
lattice displacement, U as a function of the sample temperature at different PCDs. The 
data in (a), (b) and (e) are connected by solid lines as a guide to the eye. 

Figure 5: (a) Time domain data at lowest pump fluence of 14 /zJ/cm 2 at 3K (open circles) 
along with fit to Eq. (fl|) (solid line). Inset is the FFT of the data and the fit. (a) Phonon 
amplitude, B ph vs T, (c) phonon damping parameter, ^ ph vs T (d) the phonon frequency 
v p h vs T. The solid lines are based on cubic anharmonic fit (see text) and the open stars are 
the data. 
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Figure 6: Temperature dependence of (a) u p h, (b) j p h, (c) Az/ e /_ p / l and (d) A7 e /_ p /j at 
different carrier density levels. Legends C1-C5 corresponds to 0.1, 0.2, 0.4, 0.6 and 1.4 x 
10 21 cm~ 3 . The solid lines in (a), (b) and (c) are the linear fits. The dashed lines in (d) are 
the guide to the eye. 
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